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THE  EFFECT  OF  RADIATION  FORCE  ON  SATELLITES  OF  CONVEX  SHAPE 


By  Herbert  B.  Holl 
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cooperation  and  many  helpful  discussions  of  the  mechanical  aspects  of 
the  subject.  He  also  wishes  to  thank  Mr.  James  E.  Mabry  who  carried 
out  part  of  the  numerical  computations.) 

SUMMARY 

Computations  of  satellite  position  and  attitude  will  have  to  take 
into  account  the  pressure  exerted  by  the  solar  radiation  which,  in 
general,  results  both  in  a  force  and  a  torque  acting  on  the  body. 

The  well-known  value  of  the  solar  constant  serves  for  the 
computation  of  the  integrated  radiation  force  which  is  given  the  form 

F  =  f>.  ■  A  ■  CF 

in  analogy  to  aerodynamic  forces  (A  =  projected  area;  p0  =  radiation 
pressure  in  the  vicinity  of  the  earth).  The  "radiation  force 
coefficient"  Cp  has  been  determined,  in  the  present  study,  for  several 
convex  bodies  and  was  found  within  the  limits 

0  <  CF  <  & 


Since  part  of  the  radiation  force  is  caused  by  the  reflection  of 
light  the  optical  properties  of  common  materials,  as  transparency, 
absorption,  reflectivity,  have  been  reviewed.  Their  dependence  on 
angle  of  incidence  and  wave-length  complicates  the  computation  of  Cp. 
However,  adequate  values  cannot  be  obtained  from  geometry  alone  without 
due  regard  to  the  physics  involved  in  the  reflection  process. 
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I.  INTRODUCTION 
A.  General 


Some  of  the  exterior  forces  that  will  or  may  act  on  a  satellite 
are  provided  by: 

Gravitation 

Air  Drag 

Radiation  Pressure 

Meteoric  Impacts 

Cosmic  Rays 

The  Magnetic  Field  of  the  Earth 

Orbit  calculations  are  largely  based  on  the  gravitational  and  aerodynamic 
forces.  Attitude  control  and  determination  of  position,  however,  require 
the  study  of  the  remaining  forces. 

The  object  of  the  present  report  is  to  investigate  the  influence 
of  the  radiation  pressure  of  the  sun's  light. 

When  studying  the  radiation  effect,  we  have  to  consider  separately 
the  influence  of  the  incoming  light  and  that  of  the  reflected  light.  The 
latter  depends,  in  a  large  measure,  on  the  reflectivity  of  the  opaque 
surface  material.  For  transparent  media,  the  transparency  also  has  to  be 
considered.  All  bodies,  however,  show  a  certain  degree  of  reflectivity. 
Both  the  reflectivity  and  the  transparency  are  functions  of  the  angle  of 
incidence  of  light,  of  the  wave-length,  and  of  the  refractive  index. 

These  facts  complicate  our  present  problem.  It  will  be  necessary  to 
derive  exact  formulas  for  the  force  produced  by  the  radiation.  In 
analogy  to  the  aerodynamic  equation  for  the  drag,  the  radiation  force, 

F,  acting  on  a  body  will  be  expressed  in  the  form: 

■=»  jb*  1  A  ’  y  (1) 

where  A  means  the  projected  area  of  the  body  in  the  direction  of  the 
light,  and  p0  is  the  radiation  pressure  of  the  sun's  light  in  the 
vicinity  of  the  earth.  We  can  consider  both  p0  and  A  as  given  and> 
therefore,  it  is  necessary  to  derive  formulas  for  the  coefficient,  Cp, 
which  we  refer  to  as  the  radiation  force  coefficient.  Determination  of 
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this  coefficient  Cp  is  the  first  step  in  the  studies  concerning  the 
exterior  forces  acting  on  satellites.  The  present  report  will  especially 
deal  with  the  case  of  convex-shaped  bodies,  A  second  publication  on 
concave-shaped  bodies  is  in  preparation. 

For  all  computations,  some  physical  data  have  to  be  known, 
especially  the  reflectivity.  A  few  useful  tables  are  given  in  this 
report.  For  some  general  purpose  a  set  of  tables  containing  the  values 
of  the  reflectivity  in  the  "Gaussian  plane  of  the  indices  of  refraction" 
is  planned. 

B.  The  Radiation  of  the  Sun 


In  this  report,  only  the  wave  radiation  of  the  sun  in  the 
vicinity  of  the  earth  is  considered.  Corpuscular  radiation  will  be 
excluded.  Since  the  energy  curve  of  the  sun's  wave  radiation  includes 
an  interval  stimulating  the  sensibility  of  the  human  eye,  the  sun's 
radiation  is  mainly  felt  as  light  and,  therefore,  the  expression  "light 
pressure"  is  often  used  instead  of  "radiation  pressure." 

The  total  radiation  of  the  sun  amounts  to  approximately  3.78  X 
l(r  erg/sec.  The  radiation  emitted  by  the  sun  extends  from  short  ultra 
violet  to  long  infra-red  wave-lengths.  According  to  its  absorption 
spectrum,  the  sun  is  a  GO-star  with  the  maximum  radiation  at  the  wave¬ 
length  4780  angstroms.  Diagrams  that  show  the  distribution  of  energy 
over  the  whole  spectrum  are  well  known.  (Ref.  l) 

Only  a  portion  of  the  sun's  output  reaches  the  vicinity  of  the 
earth.  The  free  space  intensity  of  the  solar  radiation  flowing  per 
unit  time  through  a  surface  of  1  cnr  placed  normally  to  the  light 
direction  of  the  earth's  mean  solar  distance  (1,496  x  10^  cm)  is  known 
as  the  solar  constant  of  radiation,  S,  and  has  the  value  (Ref.  2) 

S  =  1.97  *  0.01  cal  cm’2  min-^ 

o  i  (2) 

=  1,374  •  10&  erg  cm"^  sec 

If  the  radiation  is  to  be  considered  over  a  long  time  period,  it  is 
necessary  to  include  the  change  in  distance  between  the  earth  and  the 
sun.  Because  of  the  eccentricity,  e,  of  the  orbit  of  the  earth  the 
intensity,  J,  of  the  radiation  varies  by  7%  between  perihel  and  aphel: 
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J£ 

Ja 


(3) 


The  values  of  the  energy  of  the  sun's  radiation  at  the  extreme  distances 
are: 


Earth-Sun 

Energy  of  Sun's 
Radiation  in 

Distance  in  cm 

erg  cm"^  see-* 

Perihel 

1.4710  x  1013 

1.351  x  106 

Aphel 

1.5210  x  1013 

1.397  x  106 

The  earth-satellite  distance 

is  small  when  compared 

to  the  earth-sun 

distance.  For  example,  in  the  case  of  a  satellite  at  a  distance  of 
nearly  6  earth  radii,  the  error  in  radiation  intensity  is  only  0.12% 
using  the  earth-sun  distance  instead  of  the  sun-satellite  distance. 

In  order  to  compute  the  effect  of  the  radiation  on  a  satellite, 
the  radiation  of  subdivided  parts  of  the  spectrum  is  needed.  In  Table 
I  the  distribution  of  energy  in  the  spectrum  of  the  sun  is  given,  taken 
from  Linke's  "Meteorologisches  Taschenbuch. "  (Ref.  3)  This  table 
gives  the  sun's  energy  in  cal  cm"  2  min"l  per  100  angstroms.  The  values 
of  Table  I  are  plotted  in  Figure  1,  curve  "a."  The  scale  is 


Abscissa 

1.00  cm 

0.2  x 

10"^  cm 

Ordinate 

1.00  cm 

2-103 

O  -1 

cal  cm  J  min 

and  hence  1  cm^  of  the  graph  is  equal  to  0.04  cal  cm"^  rnin"^-.  Thus  the 
total  area  under  curve  "a"  represents  the  value  S  =  1.98  cal  cm"^  min’^. 

From  Table  I  another  graph  was  established,  Figure  2,  in  which 
curve  "a"  gives  the  amount  of  the  solar  constant  from  %  ■  0  until  a 
desired  wave-length  is  reached  or  for  any  wave-length  interval  by 
using  the  difference  in  the  values  of  two  ordinates. 

The  sun's  radiation  stems  from  a  surface  region  whose 
temperature  is  nearly  6000°  K.  As  the  sun's  radiation  sometimes  is 
likened  to  that  of  a  black  body,  it  is  of  interest  to  plot  the 
distribution  of  radiation  density  for  a  black  body  of  the  temperature 
6000°  K,  which  is  given  by  Planck's  law  (Appendix). 

JoCX  =  ij'oCX  [■*£]  (*> 
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In  Figrue  1,  curve  "b"  shows  Planck's  radiation  function  for  a  black 
body  of  6000°  K.  Curve  "b"  in  Figure  2  presents  the  portions  of  the 
total  radiation  of  a  black  body  (6000°  K)  in  the  same  way  as  does  curve 
"a"  for  the  solar  radiation. 

C.  The  Mechanical  Force  Exerted  by  Radiation 


In  the  vicinity  of  the  earth, the  sun's  light  is  to  be  considered 
as  a  system  of  plane,  finite,  electromagnetic  waves,  traveling  in  empty 
space.  By  Maxwell's  electromagnetic  theory,  an  impulse  in  the  direction 
of  the  beam  is  associated  to  each  wave.  Let  the  quantity  J  [  Watt  cm”2^ 
denote  the  intensity  of  the  light  waves  per  time  and  surface  units. 

When  the  wave  of  light  meets  the  surface  of  a  body  perpendicularly,  the 
radiation  pressure  is  given  by 

T 


p 


/C. 


(5) 


The  value  of  J  is  equal  to  the  magnitude  of  the  Poynting  vector,  E  x  H, 
where  E  and  H  are  the  electromagnetic  field  vectors.  We  measure  the 
electromagnetic  quantities  in  the  so-called  Giorgi  system,  as  employed, 
e.  g. ,  in  reference  (4).  Therefore  p0,  the  radiation  pressure  at  the 
earth's  mean  solar  distance,  is  given  by 


£  _  UljlRL  .  ^ 

/C  /C  /C 


0,  !*f"  it7 

3  '/o'* 


/*+**  *■ 

.  ±2- 

/t+K.  btC. 


(6) 


The  force  due  to  radiation  of  normal  incidence  acting  on  a  plane  surface 
element  dA0  is: 

d?  ~  a-  OC  4 0  (7) 


If  the  light  rays  form  an  angle,  Ci.  ,  with  the  surface  normal,  the  force 
is  given  by: 


oCF  ®  /v  •  o£A4  •  0(/ 


(8) 
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It  will  produce  a  translatory  motion,  which  moves  a  body  away  from  the 
sun.  If  the  material  of  the  body  possesses  a  transparency,  only  that 
part  of  the  light  which  is  not  allowed  to  pass  through  will  contribute  to 
the  radiation  force.  The  transparency  is  a  function  of  the  angle  of 
incidence  and  of  the  refractive  index  of  the  material  which  in  turn 
depends  on  the  wave-length.  It  should  therefore  be  written  as: 

T  =  T (• «■  ) 

With  a  partly  transparent  surface  material  the  force  element  is  smaller 
than  before: 


oC F  a*  jb,  •  T’fafxJ,#))  •  /t<n oo.  ^ 

In  satellites  we  encounter  shapes  like  spheres,  cylinders, 
cones,  and  planes.  Most  satellites  are  composed  of  parts  of  such 
bodies,  e.  g. ,  the  communication  satellites.  Antennas  in  use  are  of 
spherical  and  parabolic  shape.  Only  in  special  areas  the  radiation  will 
act  normal  to  the  surface.  If  the  surface  is  not  normal  to  the 
radiation,  there  are  force  components  due  to  the  reflectivity  that  give 
lateral  displacements. 

The  reflectivity  of  a  material  is  a  function  of  the  index  of 
refraction,  i.  e.,  of  the  wave-length,  and  of  the  angle  of  incidence; 
thus  we  write: 


K  -  %  ('*■  C%)/  J  . 

In  the  following  discussion  "ideal  reflection"  will  be  understood  to 
mean  that  the  beam  follows  the  ordinary  reflection  law,  and  that  it 
experiences  no  loss  in  energy  (R  =  1.) 

In  summary,  it  can  be  stated  that  the  radiation  force  acting 
on  a  satellite  is  a  function  of  the  following  quantities: 

1.  The  sun-satellite  distance,  A  , 

2.  The  size  and  the  shape  of  the  surface, 

3.  The  optical  properties:  reflective  ratio  and 
transparency  ratio, 

4.  The  angle  of  incidence  of  light,  <b  , 

5.  The  distribution  of  the  radiation  energy  in  the 
solar  spectrum. 
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Before  we  study  the  effect  of  radiation  pressure  on  selected 
bodies  of  different  shapes  we  will  give  a  survey  about  some  facts  in 
optics,  which  will  be  needed  for  a  proper  evaluation. 

D.  Optical  Properties 


a.  The  Reflecting  Power  of  Transparent  or  Dielectrical  Media 


The  oscillation  of  the  incident  light  may  be  assumed  to 
consist  of  two  equal  components  that  are  plane-polarized  in  mutually 
perpendicular  azimuths:  one  perpendicular  to  the  plane  of  incidence 
and  the  other  in  the  plane  of  incidence. 

The  following  symbols  are  used: 


E_^  ,  E/j  :  Amplitudes  of  the  electric  field 

intensity  perpendicular  and  parallel 
to  the  plane  of  incidence 


J 


Intensity  of  the  light 


i,  r,  t  :  Indices  referring  to  incident, 
reflected  and  transmitted  beam 


Angle  of  refraction 


The  amplitudes  of  the  electric  field  intensity  are  given  by  Fresnel's 
reflection  equations: 


Z  /***./$ -aft  o( 
Z  't+u-/3'&SO  q i 


(10) 
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From  these  intensities  become: 


^  ^ { <>(•& 

st+H*  Z  { 


£**■  *  (<X -ft) 
tuc  Z  (<*+/*) 

sU4+  /8o( 

l{ ' OC  ) 


^S/HC  Jot  ‘  Sp**+ 


/l***  Ytilft)  ■  A+*lCo(-{i) 


(11) 


The  reflectivity,  R,  and  the  transparency,  T,  are  defined  by: 


(12) 


(13) 


The  angle 


ft 


is  related  to  OC  by  Snell's  law: 


/A*+*'  oC 


s*  //9'C' . 


(14) 
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where  n  is  the  refractive  index. 

x 

There  is  a  special  case  for  tanofi/=  n  (Brewster's  law),  where 
becomes  zero  and  all  of  the  reflected  light  oscillates 
perpendicularly  to  the  plane  of  incidence.  For  normal  incidence  (angle 
oO  =  0,  therefore,  0=0),  the  equations  of  reflection  and  transparency 

can  be  transformed  into  the  expressions: 


R  becomes  unity  only  when  n  is  infinite,  which  is  the  case  of  an  ideal 
reflector.  Materials  with  infinite  refractive  indices  have  an  electrical 
conductivity  and  will  be  considered  in  the  next  section. 

From  the  equations  (12)  the  reflectivity  can  be  calculated. 
Figure  3  shows  the  variation  of  reflecting  power  with  the  angle  of 
incidence  for  crown  glass  (n  =  1.5):  curve  a^  for  light  polarized  in 
the  plane  of  incidence,  curve  a^  for  light  polarized  perpendicularly 
to  the  plane  of  incidence.  Since  it  was  assumed  at  the  very  beginning 
of  this  section  that,  for  incident  light,  the  two  intensity  components 
are  equal,  it  is  assumed  in  the  figure  that 


b.  The  Reflecting  Power  of  Absorbing  Materials  or 
Electrical  Conductors 


These  materials  include  metals  and  alloys.  Instead  of  the 
real  refractive  index  n  it  is  convenient  here  to  use  a  complex 
refractive  index  N  =  n-ki  =  n  (1  -  K  i) .  In  this  expression,  n  is  no 
longer  defined  in  the  same  way  as  in  Snell's  law.  This  law  is  still 
effective  when  n  is  replaced  by  N,  but  it  has  only  a  formal  meaning. 

The  angle  Q)  now  assumes  a  complex  value,  (3  =  'f'  i.  Fresnel's 

equations  for  the  intensity  of  light  change  to  (Ref.  5): 
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'Tjt  *{o(  -  e/)  Z  ^ 

\ Tj  ^z(*+9>;  y 

(17) 

f'****  l^'  ^  *(**  rf+sUHt  •tj 

JJT  6  ^***  lfap)+A^*l  lfJ{ ^(oC-PJy/sU**.  iffJ 


where  the  angles 

^ and  ^  are  evaluated  from: 

stint, 

=  -  y  -f-  f  £*  + 

✓o* 

Y- 

,0*1  f 

(18) 

l 

« 

tC  oC 

s£  st44t*?C 

j  Zo</ 

,  Z0C 

c 

The  quantities  <2  and  A  in  the  expression  for  N 

are  termed  as  follows: 

4,  =  extinction  or  absorption  index 

zjt  =  extinction  or  absorption  coefficient 


These  quantities  can  be  determined,  e.g.,  by  measuring  the  light 
intensity,  J,  transmitted  through  thin  films  of  the  opaque  material. 
The  following  formulas  are  found  to  hold: 
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with 


where 


T 

T 


X 

X 


T. 

X 


-  +rse<>c  /z, 

-  +rjtcc  /z 


-  m 
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wave-length  in  vacuo 
wave-length  inside  the  medium 
the  thickness  of  the  medium. 


(19) 


As  in  the  case  of  dielectrical  substances  the  reflectivity  is  given  by  a 
particularly  simple  expression,  if  <%  -  ^  =  0  (normal  incidence): 


n  « 


£ 


v 


fa-/)  V  'nit* 

( '*■+/)  1  y- 


(20) 


While  the  index  of  refraction  of  a  transparent  medium  varies  only 
slightly  with  the  wave-length  of  light,  the  index  for  metals  shows  marked 
dependence  on  wave-length.  A  great  number  of  determinations  have  been  made 
on  all  kinds  of  materials  and  the  data  are  published  in  many  physical 
handbooks,  e.  g. ,  (Ref.  6).  Values  for  the  reflectivity  of  light  are 
also  tabulated.  But  the  great  majority  of  the  tables  contain  only  the 
reflection  for  normal  incidence  of  light.  Values  for  oblique  incidence 
seem  to  be  very  rarely  published.  Generally,  it  will  be  noticed  from 
available  information  that  the  reflectivity  of  most  conducting  materials 
is  low  in  the  ultra-violet  and  increases  with  the  wave-length  to 
approximately  90%  or  more.  In  order  to  give  a  survey  on  the  reflectivity 
and  the  index  of  refraction,  two  graphs  were  made.  Figure  4  shows,  in 
the  Gaussian  plane,  the  index  of  refraction,  N  =  n  -  ki,  for  a  number 
of  metallic  elements;  the  values  are  taken  from  Reference  7.  The  abscissa 
denotes  the  real  part  of  the  refractive  index,  and  the  ordinate  the 
imaginary  part.  The  different  curves  refer  to  selected  metals  and  show 
how  their  index  N  changes  with  the  wave-length.  The  three-place  values 
at  the  beginning  and  end  of  the  curves  give  the  wave-length  in 
(=  10"^  cm).  In  order  not  to  overcrowd  the  graph,  only  some  samples 
of  elements  are  shown.  From  equation  20  the  reflectivity  at  zero 
angle  of  incidence  was  computed  and  plotted  in  Figure  5.  The  coordinates 
are  the  same  as  in  Figure  4  and  the  curves  give  the  reflectivity  in  0.1 
intervals.  If  the  reflection  power  with  incident  angle  other  than  0°  is 
required,  a  distinction  must  be  made  between  the  light  oscillating  in 
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the  plane  of  incidence  and  perpendicularly  to  it.  Therefore,  for  each 
angle  of  incidence,  two  diagrams  like  Figure  5  are  necessary,  where 
again  the  same  coordinate  system  as  in  Figure  4  will  be  used.  The 
refractive  index  N  depends  both  .on  the  material  and  wave-length  of  the 
incident  light.  For  each  material  and  wave-length  two  curves,  such  as 
given  for  stainless  steel  (%=  436*y<,  N  =  0.9  -  2.5  i)  and  crown  glass 
in  Figure  4  have  to  be  drawn.  From  these  curves  the  reflectivity  for 
each  angle  of  incidence  is  obtainable.  The  computation  of  the 
reflectivity  for  all  materials  would  require  much  labor  and  would  be 
hardly  economical.  However,  it  is  planned  to  work  out  diagrams  such  as 

Figure  5  for^C  =  o,  10,  20, . 80,  85,  90°,  since  this  can  be  done 

without  reference  to  specific  materials.  The  results  will  be  given  in 
a  special  report.  From  these  graphs  one  can  extract  the  values  needed 
to  obtain  the  curve  of  the  reflecting  power  for  all  possible  refractive 
indices.  Such  tables  apparently  are  not  published  as  yet,  and  it  is 
felt  that  they  would  be  of  great  help  in  engineering  work. 

Some  remarks  with  respect  to  the  refraction  indices  are 
necessary.  Physical  tables  give  the  values  of  optical  constants  with 
great  accuracy,  and  the  values  of  the  refractive  indices  of  transparent 
media  are  presented  to  seven  decimal  places.  But  the  determination  of 
optical  constants  such  as  the  index  of  refraction  requires  experiments 
which  are  among  the  most  difficult  ones  in  physics.  In  addition 
materials  are  not  easily  prepared  or  become  dull  within  a  short  time. 

Thus,  the  problem  may  become  almost  unsolvable  in  some 
instances.  In  conclusion  we  may  say  that  values  of  the  refractive  index 
with  an  accuracy  of  one  per  cent  are  desirable.  Nevertheless,  the  values 
published  for  opaque  materials  are  useful  to  get  an  idea  of  the  order 
of  magnitude  of  the  reflectivity  for  different  light  frequencies.  In 
the  case  of  newly  developed  materials  used  for  satellites,  special 
measurements  of  the  reflectivity  will  be  necessary. 

II.  THE  EFFECT  OF  RADIATION  ON  BODIES  OF  DIFFERENT  SHAPES 

Among  the  surfaces,  the  plane  surface  element  evidently  is  the  most 
important.  In  the  following  sections  we  will  consider  the  element  and 
the  plane  itself.  We  shall  also  deal  with  bodies  which  are  symmetric 
with  respect  to  the  light  direction.  Such  bodies  can  have  convex  or 
concave  form.  In  this  report  only  bodies  with  a  convex  surface  will  be 
considered.  A  future  report  will  deal  with  the  behaviour  of  the  radiation 
pressure  on  concave  bodies,  among  which  the  parabolic  and  the  spherical 
mirrors  are  the  most  important.  Finally,  cylindrically  shaped  bodies 
are  of  interest,  and  we  shall  study  here  the  right  circular  cylinder 
whose  axis  of  symmetry  is  perpendicular  to  the  direction  of  light. 


2X 
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A.  The  Surface  Element 

We  consider  a  surface  element  dA0  with  the  exterior  normal  n 
located  at  the  origin  of  a  rectangular  coordinate  system,  x,  y,  b,  with 
unit  vectors  i,  j,  k.  The  light  travels  in  the  negative  y-direction; 
after  the  reflection  on  the  surface  element,  it  has  the  direction  A  . 
From  geometric  optics  it  follows  that  the  incident  angled  ,  the  angle 
between  n  and  -j,  is  the  same  as  the  angle  between  n  and  A  ,  and  that 
the  three  vectors  are  in  one  plane,  called  the  plane  of  incidence.  This 
plane  cuts  the  plane  y  =  o  in  a  line  which  makes  an  angle,  T*",  with  the 
e-axis,  Figure  6.  V 

.  _  T^e  acting  radiation  force  we  shall  write  in  the  form 

F  =  xi  +  yj  +  e  k.  Then  the  force  element,  produced  by  the  incoming 
light  emerges  as: 


w  -/>„  ■  ({'*■)  4,  ■( /-  rc«))  j 

S  —  ^  '  /C&7  &  ■(/  -  rfrjJ  f 


The  radiation  force  produced  by  the  reflected  light  works  in  the  opposite 
direction  of  the  reflected  beam  and  becomes: 

oC  F  *  as  -8  (<*J  si  (22) 


Since 


4  sr  j  t ••  <6*3 fa  A 


(23) 


and 


/K  »  C  7*-  4*3  ^4*3^  AA4+OC  /A.  J 


(24) 


is  given  by: 
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m  - •  %  («) 


(25) 


( % t(  i  y-e^y *■«*? p /Uc+A x  Ji) # 


In  the  differentials  for  the  radiation  force,  given  by  the  equations  (21 
and  25),  the  functions  T  (#)  and  R  (et)  are  still  unknown.  Values  of 
T  (n(^),«t)and  R  (n(A),«0  for  given  N,  n,  and  c(/  can  be  calculated 
by  the  equations  (10)  through  (18).  The  characteristics  of  both 
functions  are  the  same  and  in  order  to  discuss  them  we  use  only  R.  We 
consider  a  surface  element  for  which  the  angle  interval  cia  <  d  < 
is  small  enough  so  that  the  reflectivity  is  only  a  function  of  the  wave¬ 
length  of  the  light.  The  part  of  the  reflected  intensity  which  belongs 
to  the  wave-length  %  is  written  as  R  (n(^) ,  cL  )•  where  **  refers 

to  the  intensity  of  the  incident  light.  Then  the  quantity  R  («,) 
obtained  by  the  integration  of  the  reflectivity  over  the  whole  energy 
spectrum: 


H  (*)  = 


/  "H  ['*■(*■),«]■  7X  oCX 

co 

J  7X 

0 


(26) 


gives  the  reflectivity  associated  with  the  angle  O/  in  the  above  range. 
In  order  to  master  this  integration  the  radiation  is  divided  into 
d  intervals  in  which  values  of  R  are  represented  by  R  ( #)  and 


given  by: 


Tty  U) 


Thus  R  ( ot)  becomes: 

3  «)  • 


with  /i iaQ 

4 


^ _ r  a 

Z.  3.  («)  J  ^  ** 
sh  l 


J 

(?C 


V  ** 


j  X* 


(27) 


and 
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In  the  same  way  we  introduce  a  transparency  function  T  ( eO) : 


where 


/  oCx 


6 


In  order  to  calculate  these  functions  we  can  use  tables  or  graphs,  such 
as  given  by  Figures  2  and  4. 


B.  The  Body  of  Arbitrary  Shape 

The  force  acting  on  an  arbitrary  body  is  obtained  by  integration 
of  the  forces  which  act  on  the  surface  elements.  Thus  we  have: 


/-F‘/  a  F‘ 


*  J  oC  F  6 

'  /■  4 


(30) 


Sltlfatt  4 0 

The  surface  element  dAQ  will  depend  on  two  variables,  one  of  which  we  can 
often  take  as  the  angle  of  incidence,^,  whereas  the  other  one  will  be 
determined  by  the  shape  of  the  body  we  consider.  Let  this  variable  be 
denoted  by  (  .  In  most  cases,  the  surface  element  can  then  be  written 
as:  / 


c/s  A 


aCoc,  aCj  f 


(31) 


so  that 
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For  short  we  put: 


S  («) 


Thus  the  force  takes  on  the  form: 


F*  »  p.  J  f'-  r<V]‘  <? M  •/tovol  dC-oC 


But  the  values  of  the  functions  T  (*)  and  R  («M  are  given  only  for 
angle  intervals  f  £  06  <  .  Therefore,  we  have  to 

sum  the  integrals  over  all  angles  eo  : 

»  r% 

F*-  f>.  Z.  J 

•6*/  oC 

■4”( 


where 


r.  -  - - —  /  T(«)tu. 

«*  -V,  i 


(36a) 


Analogously 


- —  / 


A  A’1  A-t 


(36b) 


The  integral  becomes  dimensionless  by  dividing  through  the  projected 
area  C*44**]  as  seen  ^rom  c^e  direction  °f  the  incoming  light: 

r-lft-j-i  /Vrj- 

*■/  < 

-4-i 

C )  ( c( )  •<t4?oC  dC(X/ 
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In  the  same  way,  the  radiation  force  of  the  reflected  light  becomes: 

f*-.  s<*)- 

tM/  (38: 

•  (  't+H.'p/iud.&vi  2  «(/&,)  eCo(/. 


We  introduce  the  notation  Cp  and  write  the  force  components  as: 


*  A  M 


a  -  ^  c; 


The  coefficients  Cp  are  given  by: 

)> 


y  ^ 

—  Z  j  (l-T^)-<j(«)-  *»*<**' 


•id  *L 
-6-/ 

/  ~  l  \  l 


c;  -  \i  te)'  +  (c:r+(ctf 


V  *0, 

-r  Z  J  4  °*«  *"t**u4*L< 

>«  V/ 

X> 

_L  2  ]  6  w  •  \  •  *r,<*  et,tL 

A  h>l 

» 

— 2  I  4  Otop/Huioceix 
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The  equations  (39),  F  =  Cp  x  A  x  pQ,  (41) 

have  the  shape  of  the  aerodynamic  equation  for  the  drag: 

D  =  Cq  x  A  x  q,  (42) 

(Cq  =  drag  coefficient,  q  =  dynamic  pressure). 

Therefore,  cp  might  be  termed  the  "radiation  force  coefficient";  it  is 
a  function  of  the  shape  and  of  the  optical  properties  of  the  body 
considered. 

C.  The  Surface  Element  of  a  Symmetrical  Body 

In  the  following  we  consider  only  such  bodies  which  are 
symmetrical  relative  to  the  light  direction,  in  the  sense  that  the 
individual  force  components  normal  to  the  incident  light  cancel  each 
other,  e.  g. ,  spherical  and  parabolic  surfaces  and  the  cone  and  cylinder 
in  suitable  orientations.  The  only  force  component  is  then  in  (-j) 
direction.  One  exception  will  be  the  plane  surface  inclined  towards 
the  light  direction. 

Let  dA  in  Figure  7  be  a  surface  element  of  any  symmetrical 
body;  let  the  p¥ane  of  incidence  be  the  drawing  plane.  The  element  of 
the  projected  area,  dA,  becomes: 

OCA  ~  oC /)j  ' /***  rts 

Then  the  acting  force,  dF,  is  given  by: 

eC  F  =  oCPy  +  d.  Fy  ^  fa  .  at  /)  •  C  p 

-  fa  o0Ao  -tn*  + f>9  am Zo( 

i  r 

The  radiation  force  coefficients,  Cp,  and  Cp,  are: 

« 

Cp  »  /  -  T(*) 

Cp  =  («)  mjt+9  £  • 


(44) 
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Special  values  for  the  radiation  force  coefficient  are  determined  as 
follows: 


Opaque  Materials 

1.  Medium  with  total  absorption 

ft  (tl)  s  0  f  T  («)  •  O 

/ 

2.  Medium  with  ideal  reflection 

RC«)  =  I  T(«)  *  o 

Cp  -  t  +  c»li<n  a  Z  /**>*<>(.■ 

3.  Medium  with  partial  reflection 

o  <  <i ,  t («)  *  O 

C„  a  1+  R(ct)  •  . 

h 

Transparent  Materials 

4.  Medium  without  absorption 

XM  +  Tfr)  -  / 

(*  s  l-T(«)  +  TIC*) 

s  TlC«)-(  / 

-  z  H,  («)  ■  ■ 

From  the  above  four  cases  it  can  be  seen  that  the  radiation  force 
coefficient  attains  the  value  2  only  at  normal  incidence  of  light  and 
for  ideally  reflecting  materials  (case  2).  At  an  incident  anglefl^  =  45° 
the  reflected  light  yields  no  force  component  in  the  direction  of  the 
incident  light.  We  then  have  Cp  =  1  for  all  opaque  materials.  For 
angles  of  incidence  pf  ?  45°  all  reflecting  materials  have  a  force 
component  opposite  to  the  light  direction;  Cp  is  smaller  than  unity  and 
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goes  to  zero  With  GO  +  —  . 

D.  The  Plane  Surface  (Solar  Sail) 

We  consider  a  plane  disc  which  may  serve  as  an  idealized  model 
of  a  solar  sail,  the  flight  dynamics  of  which  will  be  discussed  briefly. 
The  position  of  this  plane  body  is  indicated  in  Figure  8. 

In  the  absence  of  reflectivity,  the  disc  will  be  accelerated 
in  the  direction  of  the  sun's  light  independent  of  inclination. 

Therefore,  only  drag  is  present.  A  reflecting  disc  also  has  a  lift 
component.  From  the  equations  (21)  and  (25)  where  y  can  be  taken  as  zero 
the  acting  force  on  a  solar  sail  becomes: 

F  «  Fy  +  Pz  s  A  C0  1  *)  («) 


with 


*  /  -  T(«) 

s  %(*)•*+**« 

a  %(«)•  '*+“'*’*  * 


(46) 


The  course  of  the  solar  sail  is  determined  by  the  lift/drag  ratio  which 
gives  the  course  angle  : 


For  opaque  materials 


oC 

An 


cj_  _  <»> 

c;  "  /  +  n(<)c*a< 
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For  transparent  materials  without  absorption 


Z2**  OC  (48) 


For  an  opaque  plane  sail,  the  values  of  the  radiation  force  coefficients 
^  and  Cq  as  functions  of  R  (e^)  are  given  in  Figure  9.  The  maximum  of 
the  lift  Fb  is  found  by  differentiation  with  respect  to  tU  of  the 
expression'. 

F£  a  -  A 

while  R  is  kept  constant.  The  maximum  of  Fg  occurs  as  the  angle 

oC  =  3S.234 

without  dependence  on  the  value  of  the  reflectivity. 

It  is  worth  mentioning  that,  in  the  case  of  a  transparent 
plane-parallel  plate,  we  get  a  torque  due  to  the  parallel  displacement  of 
the  light  beams. 

E.  The  Spherical  Surface 

We  consider  a  spherical  satellite,  the  right  hemisphere  of 
which  is  illuminated  by  the  sun's  light  (Figure  10).  The  radiation 
incident  on  an  arbitrary  annular  section  dA0  is  uniform  over  the 
circumference  of  the  section.  We  get  this  ring-shaped  element  when  we 
consider  a  surface  element  da0  which  is  given  by: 

oCa,t  **  A*cUl  cC  jCk . 

By  integration  we  have: 

»r  ir 

s  J  oCq,0  =*  J  FZsi*U'ot cC*d>/b 

(1*0  ft:*  (49) 


7YO  *■  R(U) 
"Uac  OtG  »  • 


et  oL 


/  t4,  cmRoi 

<<c  *  oO  . 


Z^T" 'oli-ct  cid  5  G  (*t)  •  cCct 
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The  radiation  force  acting  on  a  spherical  calotte  with  the  half-opening 
angle  eCy  is  obtained  by  using  the  equations  (35)  through  (40): 


F 


0 

Fri*Fr’  ‘  A 


V>  of. 


V  f  _  _  /  v 


Fy*-*  Z  2  I*1' 

JLaf  *A-f 

i  r 

The  radiation  force  coefficients  Cp^  and  CpY  are  found  as 

»  t  ** 

C*’  a  ——  ■  J  ( l-T.) 

*v  /4««  .  a/ 

7  V  4*  i  \m( 

v  «r. 


-  * 

_£ _  J  ^Au+tCcnZKonoCoiof/ 


e£*/  (✓ 

y  i«/  V/ 


Here  the  functions  Tfc  and  are  given  by  the  equations  (26)  through 
(29),  and  (36a,  b) .  In  the  following,  we  consider  opaque  materials  and 
then  have  Tt  =  0.  Of  special  interest  are  materials: 

1.  With  total  absorption,  R  =  0 

2.  With  ideal  reflection,  R  =  1.  We  obtain  for  the 
acting  forces: 


/'  F  X  Fy 


■  A  •  Ci 


f  9 

s  f°e  ( •  C P 


(  C  -  — - -  oCcC 

c  .  •  i,  J 

ry  0(y  0 
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*•  F  -  'y 

y  =  />.  /)■  c£  =  /> .  (***.*,) *r  c* 

y 


1  1 

w 

C*9  MU++  V  ~J 

PY 

loYy  ' 
p 

*  I 

(54) 

l0(y  L 

1 

'  / 

v  -  A  ) 

-  /tav  t(y  -  o(y  +  /  J 

X  /<S #9 

F  » 

A  ^ 

•  p,  A  (/ +  <***■«„)  (55) 

In  Figure  11  the  values  for  the  radiation  force  coefficient  in  the  above 
two  cases  are  plotted  in  their  dependency  on  the  half -opening  angle  cO#. 
In  the  case  of  a  hemisphere  with  ideal  reflection,  c£  becomes  zero  and 
the  effect  of  the  radiation  is  equal  to  that  of  a  hemisphere  with  total 
absorption.  But  this  result  is  true  only  with  the  assumption  Rr  =  1. 

In  the  interval  0°  £  tfnp  £  45°,  spherical  surfaces  with  a  reflectivity 
0  <  R  <  1  will  be  characterized  by  curves  for  the  radiation  force 
coefficient  between  curves  1  and  2  of  Figure  11.  Values  for  CFy  for 
spherical  surfaces  with  larger  opening  angle,  ~±L  <  of  ^  —  ,  are  not 
limited  by  these  two  lines.  ^  & 


The  values  for  the  reflectivity,  R  (n  (%),&),  are  not 
available  today.  Therefore,  we  are  unable  to  calculate  the  force  acting 
on  a  metallic  sphere.  But  we  can  calculate  the  limits  for  all  kinds  of 
materials.  The  upper  limit  will  be  obtained  with  a  surface  material 
excluding  any  force  components  opposite  to  the  light  direction.  The 
lowest  values  of  (cO)  appear  when  the  component  of  the  force  Fr  in 
direction  of  the  incident  light  is  zero  and  when,  in  addition,  the 
component  opposite  to  it  is  acting  in  its  full  amount.  The  limiting 
values  of  Cpy  are  given  by  the  curves  3  and  4  in  Figure  11.  On  the  body 


24 


which  is  represented  by  curve  3,  the  area  of  the  surface  between  the 
angles  of  incidence  cO-  0°  and  <&  =  45°  consists  of  an  ideally  reflecting 
material  while  the  zone  from  CO  =  45°  to  CC  =  90°  is  assumed  as  totally 
absorbing.  We  get  the  minimum  force  if  the  inner  part  of  the  sphere's 
surface  has  total  absorption  and  the  outer  zone  possesses  ideal 
reflection.  When  the  entire  hemisphere  is  illuminated,  the  values  of 
Cp^  are  limited  by: 

0.75  <  C_  i  1.25 
—  Fy 


The  Cylindrical  Surface 


We  consider  1)  a  right  circular  cylinder,  with  its  axis 
of  symmetry  in  the  a-direction,  as  indicated  in  Figure  12a, and  2)  a 
cylindrical  surface,  which  might  be  part  of  an  arbitrary  body  (Figure 
12b).  For  symmetry  reasons  the  resultant  of  the  force  components  in 
x-  and  z -directions  becomes  zero.  The  y-component  can  be  calculated 
from  the  equations  (21)  and  (25).  The  anglg  oClies  in  the  x-y-plane. 

In  the  case  of  Figure  12a  we  have  »  -j*  .  The  surface  element  dA0 

is: 


ft  '/t  *  &Cc C 


If  Tt  =  0,  the  force  of  the  incoming  light  becomes: 

Fy  St  jt>4  •  J  A+V  U  doc 

ss,  2/i£/U*+o(y  •  •  — —  J /tjy  o(  aCo(/ 

f*  ”  9 

cl  =  }*»*»**  s  t 

7  ***  , 


The  reflected  light  contributes  the  force: 


F ~  * 


■.  z  l  ■*. 


'‘h 


=  S  • 
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In  Figure  13,  curves  1  and  2  give  CFy  for  Rr  =  0  and  Rr  =  1.  With  the 
same  arrangement  as  used  with  the  spherical  surface,  we  get  the  limits 
for  a  cylindrical  body.  Curve  3  in  Figure  13  represents  the  case  where 
Rr  =  0  from  0^  *  0°  to  CO  =  45°  and  where  Rf  =  1  from  oO-  45°  to  cC  =  90°. 
Curve  4  shows  the  opposite  case  with  respect  to  the  reflection.  For  the 
whole  cylinder  we  get  the  limits: 

0.862  CFy  4  1.471. 

G.  The  Conical  Surface 


Here  we  deal  with  two  cases  only:  (a)  where  the  geometrical 
axis  is  perpendicular  to  the  sun's  radiation, and  (b)  where  the  cone's 
vertex  points  towards  the  sun.  (Figures  14a  and  b.) 

a.  Cone  Axis  Perpendicular  to  the  Sun's  Radiation 

The  surface  element  dA0  at  point  P  is  a  rectangle  formed 
by  the  slant  height  ds  and  by  r.  djf  ,  where  r  is  the  radius  of  the  circle 
through  the  point  P.  (Figure  14c.)  The  slant  height  is  counted  from 
the  vertex;  thus  we  have  r  =  s  •  sinAJ,  and  the  surface  element  dA0 
becomes: 

at  **  *  A' A+tcU dii -dX  (61) 

From  Figure  14a  we  calculate  the  exterior  normal  n: 

nH  SB  CJ  /HU  ’t/UfvCJ  C&X T  6  (62) 

for  which  we  had  previously  obtained  the  general  relation: 


m,  a  +, 4*7*  j  +/**?£’ * \XuoC  Xi 


(24) 
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By  comparison  of  the  coefficients  we  get: 


Si/*** 


/t<n  d 


/t*v 


CO  •  s€++€> 


st*+cOt 
sd>ut€/  CO 
sisU*  oC 

/€^9  6)  ■  y. 


In  the  case  of  spherical  surfaces  as  well  as  for  hyperbolic  and 
parabolic  forms,  it  is  practical  to  integrate  the  force  produced  by 
the  radiation  over  the  incident  angle  of  the  light.  Here  the  third  of 
the  equations  (63)  shows  us  that  Oi/  is  a  function  of  the  half  opening 
angle  CO  and  of  the  angle  Jlf  .  When  a  conical  surface  is  only  a  part 
of  a  circular  cone  as  in  Figure  14,  the  size  and  the  shape  are  determined 
by  a  suitably  chosen  value  of  the  angle  .  Now  we  insert  the  results 
(63)  into  the  equations  (21)  and  (25).  The  differential  of  the  radiation 
force  for  the  incoming  light  is  then  found  as: 

i 

oiF*  ss  ~  4S9  jf- Sisk  cCXj  f  (64) 

where  the  angle  jf  plays  the  role  formerly  assigned  to  the  angle  0(/. 
Similarly,  for  the  reflected  light: 


oC  s  -pa  'i'+W  %(«) 6*?X‘St*Ci atX  • 

(Z  *nlco  /*<*>  *Jf  -  /)  j  +• 

z,  /X44+CC  4*94?  449  X  £ 
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Formal  integration  of  the  differential  equation  64  yields: 

.  *  fX+ 

F  M  -  S1A4+CO  CJ9CJ  Jsi  Zj 

sim*  *•/  X 


yi 


+1 


»  ,jg. 


(66) 


s  -  fa ' -j—  '/U4*6)C*9a  £  i  J  ChTJ.cvXdX  J 

tmf  «£/ 


with 


X 


s  0 


and 


/  mu, (£**.) 
J  +  |  /*+?&>  j . 


where  oC  is  in  the  interval 


+  \  *>  /  / 

—  >  • 
b  ==  *  - 


From  Figure  14d,  we  can  see  that  the  projected  area  of  a  cone  is 
si*1, si*** 6/ e+o  CJ  »  that  of  the  cone  segment  considered  then  becomes: 

A  m  s*** jCy  • 

Putting  the  last  expression  into  equation  (66^ 

F*  -  f  o-V-XdXf  (67) 

^  ^  "iml 

From  equation  65,  it  can  be  seen  that  force  elements  acting  on  surface 
elements,  which  are  symmetric  with  respect  to  the  y-z-plane,  cancel  each 
other,  so  that  the  i-component  of  the  integrated  force  in  zero.  As  a 
consequence: 

’  “•  -L 


be  i  r*. 


C<  4  A*9  4*9  9X  -  ?  * 

"*  Z>  '****’£*)  A+9C4  *X  J  • 


(68) 
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From  these  equations,  we  get  the  following  expressions  for  the  radiation 
force  coefficients: 

/  J?  r 


X  ( ['***•?] 


^  Z  X. 


V  r  IW 

C*  =  — r-— r '  5  /  T  t******"^ (§*y+2h 

cr  _^xl 

V  r  ,  X/t-t  «Sf, 

SS  ^  l  V  *5  \'**t*z£  X  I 

X  ’  i  l 

JA-( 


With  an  ideally  reflecting  full  cone  the  coefficients  become: 


C  «  / 
7 


X  a 

7 


it  -  / 

j 


r4  a  I.  ~ 

h  2- 


/i/nvU  •  /t&  Co 


The  total  radiation  force  then  is: 


F  =  -p.A(^f —  j  +  ■?• 


3X 
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These  forces  move  the  cone.obliguely  away  from  the  light  source.  The 
course  angle  06 L  in  the  (  -j ,  -  k)  -  plane  is  given  by 


"tZ*c  X-  a  - -  7 CO 

c  9 


(72) 


b.  Cone's  Axis  Pointing  Toward  the  Sun 


The  radiation  illuminates  the  lateral  area  AQ  (Figure  14b) 


A*  ■=  A  •  T'/i  s 


AgJT 


,  (73) 

The  half-opening  angle  CO  and  the  angle  of  incidence  oC  are  connected  by: 

oC  S  CO 

If  Tt  =  0,  the  force  of  the  incoming  radiation  becomes 

f*  -  -/*  —  '  **«}  . 

1  0  a  r°  ~ 


(74) 


where  , 

cf 

r 

For  symmetry  reasons,  the  reflected  light  has  only  a  force  component 
in  the  j -direction,  and  we  have: 

=*  -  ps  A1  T'  Tifa)  a+pA* 


(75) 


with 


/  +  c&ii  0(  =  a  ZeC  * 


r 


In  the  case  of  ideal  reflection,  R  (c£)  =  1, 

C  r~  ~  l  -f  *09  Zd  ®  A  , 


(76) 


0  < 


<  z 


The  limits  are: 
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Long  tapers,  due  to  the  reflection  force  component  opposite  to  the 
light  beam,  have  a  small  resultant  force  component  in  light  direction. 
As  an  example,  for  a  slender  taper  with  s/r  =  5,  Cif  becomes  &  78.5°, 
and  Cfy  &  0.08,  if  R  ^  1  (highly  reflecting  Surface).  When  the  cone 
angle  oC  is  larger  than  45°,  Cpy  becomes  smaller  than  unity. 

H.  The  Parabolic  Surface 


We  consider  a  parabolic  body  as  sketched  in  Figure  15a.  As  in 
the  case  of  the  sphere,  we  select  an  arbitrary  annular  section  dAo  and 
a  surface  element  da0  on  it.  Let  p  be  the  ha If -parameter  of  the  parabola. 
The  shape  then  is  given  by  p,  and  the  size  is  determined  by  the  half- 
opening  angle *0  ^  , 

For  the  derivation  of  the  pertinent  formulas,  we  consider  Figure 
15b.  The  equation  of  the  parabola  has  the  form: 


=  -Z 


ry 


and  as 


P 


(X/ 


(77) 


(78) 


the  radius  of  curvature  of  dA0  becomes: 


/;/  -  / 


D 


VjL 


/ 


/Cits  *o< 


/ 


(79) 


Thus  the  surface  element 


oCa,e  r  pact  •tcLfi  st  pz>  — aCccacft 


(80) 


The  ring-shaped  element  dA0  is  then  obtained  as: 


^c< 


(81) 


We  exclude  the  infinite  parabola  so  that  OC  remains  smaller  than  , 

z 
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By  using  the  equations  (35)  through  (40),  the  radiation  forces  are  found 


F  ■  ‘r'rr-  £ 

J  r  ,  7  ** 

=  o-rt)[-  ~^r] 


A"* 

s  Z  J  3 


C+oZaC 


1st  <>., 
V  «/L 


*&>  V 


-  Z(>ir^l  J  ?!£&)** 

As/  "fi./ 

*  t>irp.Z  %*  [-*■£>**«  +’zrr]  • 

Am/  k  *•  “  *V  , 

A  “I 

The  projected  area  A  of  the  parabolic  body  is  given  by: 

/4  =  j»lF  =  />  *z*»  V*,  •  7~ 

Thus  we  have  for  the  radiation  force  coefficients: 

C  p  a  —  ■  •  2  —  7 

7  77/  i  <**»*•<  L 

We  determine  C„  for  different  cases  of  R  and  T:  «•/ 

Fy 

Opaque  Material 

1.  Material  with  total  absorption 

T 9  o  z *0 
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2.  Medium  with  ideal  reflection: 

7"  3  0  ,  «  / 


y  <*v 


3.  Medium  with  constant  reflectivity: 

r  =  o  ,  o  <'R  <  t 


—  -  A*r**t*v  I 


(87) 


Cr*  S  — - -  /-  H"  R*  /lu*c#)0(  •+ 

y  l 


Tt-  (  -I^A  (88) 


/4#iV 


J 


M-/ 


Transparent  Material 

4.  Medium  without  absorption 
R  +  T  =  1 

C  =  * 


Fy  -u^^v  k„ 


J  'ft/l  !--&<■  cnu  I 

1  Joi 


1*1 


(89) 


i-/ 


From  equation  (87)  it  can  be  seen  that  in  the  present  case  is 
bounded  by:  ^ 


0  <  CF  <  ^  . 

* 


r 


III.  CONCLUDING  REMARKS 


These  studies  were  made  for  application  to  computations  of  both  a 
satellite's  attitude  and  trajectory,  which  are  to  take  into  account 
the  effects  of  the  radiation  pressure. 

The  radiation  force  acting  on  a  body  is  a  function  of  certain 
physical  properties,  of  which  the  reflectitivy  of  the  body's  surface  is 
the  most  important.  The  force  F  can  be  expressed  in  the  form: 

F  =  p0  x  A  x  Cp 

"A"  means  the  projected  area  of  the  body  and  p0  is  the  radiation 
pressure  in  the  vicinity  of  the  earth.  The  function  Cp  was  introduced 
in  analogy  to  the  aerodynamic  expression  for  the  drag,  and  may  be 
called  the  "radiation  force  coefficient."  The  purpose  of  this  report 
was  to  derive  accurate  formulas  and  limits  for  the  radiation  force 
coefficient  for  certain  convex  bodies.  For  finite  bodies,  the  result 
is  that  Cp  is  limited  by: 

0<CF  =  2 


Plane  Surfaces 

0<Cp  *  2 

Spheres 

0.75  =  Cp  ^  1.25 

Cylinders 

(Illuminated  perpendicularly  to  the  axis  of  symmetry.) 

0.862  f  CF  ^  1.471 

Cones  in  Special  Positions 
(See  Figures  14a  and  14b) 


0  <  CF  <  2 


Parabolic  Bodies 


34 


IV.  APPENDIX 

The  total-radiation,  E,  of  a  black  body,  coming  from  1  cm^  of  the 
surface,  is  given  by  Stefan-Boltzmann's  law,  which  states  that  the 
total  rate  of  radiation  varies  as  the  fourth  power  of  the  absolute 
temperature  T,  that  is: 

E  -e-T4. 

T  is  the  absolute  temperature,  and  a  constant. 

The  rate  of  emission,  E  ^  ,  of  the  radiation  of  any  particular 
wave-length^,  which  passes  in  1  second  through  each  cross  section  of 
a  radiation  cone  of  the  aperture!!  =  1,  is  given  by  Planck's  equation: 

=  /C,  X'r  ( e.  *X/zT_/)~/ *CZ 

Here  the  notations  are: 

x: ,  = 


c  =  the  velocity  of  light  in  vacuo 

k  =  the  gas  constant  for  1  molecule  (Boltzmann's  constant) 
h  =  Planck's  quantum  of  action;  with  the  values: 
c  =  299  850  km  sec"l 
k  =  1.372  Watt  sec  degree"^ 
h  =  0.655  x  10'33  Watt  sec3 
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:1  = 


5.9  x  10*6  erg  cn)2  sec-l  B  5,9  x  IQ’13  Watt  cm2 
=  1.42  x  10"^3  cal  cm2  sec*^ 


C2  =  1.43  cm  degree 


Q  = 


5.8  x  10'D  - - ££S - r 

cm*  sec  degree 


5.8  x  10'12  Watt 

cmz  degree4 


1.39  x  10*12 


cal 


cm*  sec  degree^ 


The  connection  between  Planck's  and  Stefan-Boltzmann's  equation  is 
given  by  integration  of  Planck's  equation  over  the  whole  spectrum,  as 
follows: 


/  e *  -  /c  /  *1, - -  di  =  —  r* 

This  is  the  energy  within  the  solid  angleJl  =  1.  (Reference  8). 
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TABLE  I 

THE  DISTRIBUTION  OF  ENERGY  IN  THE  SPECTRUM  OF  THE  SUN 
Extraterrestical  Radiation  Intensity  of  the  Sun  i0^  [in  10“3  cal  cm"2 
min'lj taken  from  M.  Nicolet,  Archiv  Meteor. ,  Geophys.,  Bioklim.  B  3 
(1951)  209  for^A*  =  100  AE.  Values  in  Smithsonian  Scale  1913  minus 
2.4%,  Solar  Constant  =  1.98  cal  cm_2  min_l. 


x[j*l 

°ox 

%[/*] 

1 — j 

°8X 

%[f-l 

0  ox 

0.20 

- 

0.50 

29.7 

0.80 

16.5 

1.50 

4.3 

0.21 

- 

0.51 

29.7 

0.81 

16.3 

1.55 

3.9 

0.22 

0.2 

0.52 

27.7 

0.82 

15.9 

1.60 

3.7 

0.23 

0.5 

0.53 

28.9 

0.83 

15.6 

1.65 

3.2 

0.24 

0.6 

0.54 

28.9 

0.84 

15.3 

1.70 

3.0 

0.25 

0.8 

0.55 

27.6 

0.85 

15.0 

1.75 

2.7 

0.26 

1.9 

0.56 

28.8 

0.86 

14.7 

1.80 

2.4 

0.27 

2.4 

0.57 

28.4 

0.87 

14.4 

1.85 

2.3 

0.28 

2.2 

0.58 

28.4 

0.88 

14.1 

1.90 

2.0 

0.29 

5.1 

0.59 

27.5 

0.89 

13.8 

1.95 

1.9 

0.30 

5.9 

0.60 

27.2 

0.90 

13.5 

2.00 

1.7 

0.31 

8.8 

0.61 

26.3 

0.91 

13.2 

2.05 

1.6 

0.32 

10.3 

0.62 

25.7 

0.92 

13.0 

2.10 

1.4 

0.33 

12.6 

0.63 

25.1 

0.93 

12.7 

2.15 

1.3 

0.34 

12.4 

0.64 

24.5 

0.94 

12.4 

2.20 

1.2 

0.35 

13.1 

0.65 

24.1 

0.95 

12.2 

2.25 

1.1 

0.36 

13.5 

0.66 

23.6 

0.96 

11.9 

2.30 

1.0 

0.37 

13.5 

0.67 

23.0 

0.97 

11.7 

2.35 

0.9 

0.38 

13.3 

0.68 

22.5 

0.98 

11.4 

2.40 

0.8 

0.39 

14.2 

0.69 

21.9 

0.99 

11.2 

2.45 

0.7 

0.40 

22.1 

0.70 

21.2 

■B9 

11.0 

2.50 

0.7 

0.41 

25.1 

0.71 

20.7 

■HI; 

10.0 

3.00 

0.37 

0.42 

25.2 

0.72 

20.2 

1.10 

8.7 

0.14 

0.43 

23.5 

0.73 

19.8 

1.15 

7.8 

mm. 

0.44 

27.4 

0.74 

19.0 

1.20 

7.2 

■ 

0.03 

in 

30.0 

0.75 

18.7 

1.25 

6.6 

0.01 

m 

30.9 

0.76 

18.4 

1.30 

6.1 

0.47 

30.8 

0.77 

17.8 

1.35 

5.6 

0.48 

31.2 

0.78 

17.4 

1.40 

5.1 

0.49 

29.0 

0.79 

16.9 

1.45 

4.7 
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Figure  6 


The  Surface  Element  dA0 
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Figure  7 


The  Radiation  Force  Acting  on  a  Surface 
Element  dA0  of  a  Symmetrical  Body 


The  Flight  Characteristic  of  a  Solar  Sail 


The  Radiation  Force  Acting  on  a  Spherical  Body 


The  Surface  Element  da0  of  a  Spherical  Eody 


Figure  12a 
Figure  12b 


The  Radiation  Force  Acting  on  a  Cylindrical  Body 
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Figure  14a 


The  Radiation  Force  Acting 
on  a  Right  Circular  Cone 


[Figure  14b  The  Radiation  Force  Acting  on  a  Right  Circular  Cone 
Figure  14c  The  Surface  Element 


ojected  Area  of  a  Right  Circular  Cone 


Figure  15& 


The  Surface  Element  da0 
of  a  Paraboloid 
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